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Abstract 

We give a constructive proof for the existence of an A-dimensional Bloch basis which is 
both smooth (real analytic) and periodic with respect to its d-dimensional quasi-momenta, 
when 1 < d < 2 and N > 1. The constructed Bloch basis is conjugation symmetric when 
the underlying projection has this symmetry, hence the corresponding exponentially localized 
composite Wannier functions are real. In the second part of the paper we show that by 
adding a weak, globally bounded but not necessarily constant magnetic field, the existence of 
a localized basis is preserved. 
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1 Introduction and main results 


1.1 Generalities 

The Wannier functions (as bases of localized functions spanning subspaces corresponding to energy 
bands in periodic solids) have been playing a central role since their introduction in 1937 in both 
qualitative and quantitative aspects of the one-electron theory of solid state physics. In particular, 
they are the key ingredient in obtaining effective tight-binding Hamiltonians, a topic related to the 
Peierls-Onsager substitution. After the seminal paper of Marzari and Vanderbilt [MV) , they have 
become a powerful tool in ab initio computational studies of electronic properties of materials (see 
[N2] . [MMYSV] . [FMP] and references given there). 

The existence and construction of exponentially localized Wannier functions (ELWF) has been 
one of the fundamental problems in solid state physics. Let us remind the reader that prior to 
the crucial observation made by Thouless |Th) that the very existence of ELWF implies that the 
quantum Hall current of the corresponding band vanishes, the fact that ELWF might not exist 
was largely overlooked by physicists. 

It turns out that due to several mathematical subtleties, the proof of exponential localization 
is difficult and depends on d (the dimension of the configuration space) among other things. 
Accordingly, the progress in proving the existence of ELWF in all cases of physical interest was 
slow. For example, the existence of ELWF for composite energy bands of time-reversal invariant 
Hamiltonians in two and three dimensions has been only recently achieved [P], [HPCMM] . 

The first rigorous result about the existence of ELWF in one dimension was obtained by Kohn 
in his 1959 classic paper m- Using a technique based on ordinary differential equations, Kohn 
considered simple bands in crystals with a center of inversion and obtained a complete solution: 
he constructed Wannier functions which were real, symmetric or antisymmetric with respect to an 
appropriate reflection and with an optimal exponential decay. Kohn’s method does not generalize 
to higher dimensions and the next major step was made by des Cloizeaux. In a couple of basic 
papers [dCllj . [dC12) in 1964 which are still worth reading, he showed that the existence of ELWF 
for composite bands in crystals of arbitrary dimension can be achieved in two steps: 

I . To each isolated energy band one associates a family P(k),k G of rank N orthogonal 
projections {N = 1 for simple bands and 1 < iV < oo for composite bands) which are Z'^- 
periodic and jointly analytic in a complex tubular neighborhood / of M'^. These projections 
live in a separable Hilbert space T-L. 

H . Show that the vector bundle B := {(k, v) : k G v G RanP(k)} is trivial (= the 
subspace admits a continuous orthonormal basis). More precisely, one has to show there exist 
N vectors {Sj(k)}j^]^ C Ji which are continuous, Z'^-periodic and form an orthonormal basis 
of Ran(P(k)) for all k G R'^. 

Step 1. is quite general and it is based on the Bloch-Floquet-Gelfand transform and analytic 
perturbation theory. It works for all standard Schrodinger operators with periodic potentials, 
including their discrete and pseudo-differential variants (see e.g. |dCH) . [N2] . [FMP] . |dN-G] as well 
as section [L2] below). 

Step H. is a subtle and more demanding mathematical problem, due to the fact that vector 
bundles can be non-trivial. A well known example is the tangent space to any even-dimensional 
unit sphere: the existence of a smooth orthonormal basis would contradict the Hairy-Ball Theorem, 
des Cloizeaux realized that in proving step H the obstructions have topological origin and showed 
that they are not present for simple bands {N = I) in the absence of a magnetic field or when the 
periodic electric potential has a center of inversion. The restriction to crystals with a center of 
inversion was removed in 1983 by Nenciu [NO] who proved by operator analytic methods that at 
the abstract level, for arbitrary d and TV = 1 the obstructions to the triviality of the vector bundle 
B are absent if P(k) satisfies 

eP(k)e = P(-k) (1.1) 
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for some antiunitary involution 9. The proof is constructive and leads to optimally localized, real 
Wannier functions |N2) . A simpler proof was later given by Helffer and Sjostrand [HSj : however, 
it seems that it does not provide the optimal exponential decay for the Wannier functions. 

For a Schrddinger operator with a real symbol (i.e. without a magnetic potential), (II.ID 
is insured by the fact that it commutes with complex conjugation i.e. it obeys time reversal 
symmetry. Concerning the case > 1, it has been sugested [N2] . |NN2) that exploiting (11.111 in 
the context of characteristic classes theory in combination with some deep results in the theory 
of analytic functions of several complex variables might lead to the existence of ELWF. Indeed, if 
d = 2, 3 and if P(k) satisfies (11.11) . the triviality of the bundle was proved by Panati [P], providing 
a non-constructive proof of the existence of ELWF for all cases of physical interest |BPCMM] . In 
higher dimensions, examples of non-trivial bundles have been given in |dN-G] . 

While concerning the non-constructive existence of ELWF the situation is satisfactory, there 
are still many interesting questions left open in the area. Below we list a few of them: 

1. As it stands in [P], the triviality of the bundle generated by P(k) for d = 2,3 and iV > I 
is an abstract existence result, but at least for computational aspects, a constructive proof 
is needed. An important step in this direction has been very recently made in |FMPj : the 
authors construct real valued Wannier functions which decay faster than any polynomial. 

2. All the results outlined above (as well as the very definition of Wannier functions) require 
strict periodicity of the underlying Hamiltonian. However, as stressed by many authors 
(see e.g. |KOj . |GK) . [Ki] ) results on the construction of exponentially localized bases in 
non-periodic (or at least nearly periodic) systems are highly desirable. 

3. For theoretical as well as computational purposes it is highly desirable to construct local¬ 
ized Wannier functions inheriting as much as possible the symmetries of the underlying 
Hamiltonian, e.g. time reversal invariance should lead to real valued Wannier functions. 

4. At least from a computational point of view it is important to construct Wannier functions 
with best localization properties. In this context, one has to define precise localization 
criteria and to discuss the relations between them [MVl IGNNl IMPl IPP) . 

5. In dimension one, Wannier functions can be constructed as the eigenvectors of the position 
operator restricted to the range of the Riesz projection of the isolated band |NN1| . Gan this 
method be extended to higher dimensions? Some recent progress on the so-called ‘radial 
localization problem’ has been made by Prodan [Prj . 

The aim of our paper is to add some results concerning the first three questions above. Our 
first main result is a construction (assuming (II.Ill I of real valued ELWF for d = 2 and N > 1. Our 
method is quite different from the one in [FMPj : in particular, it does not require a preliminary 
reduction to N = 2, and builds on the operator analytic approach in [NO) . |N2) . Note that the 
results of our paper are not applicable to the so-called ‘fermionic’ case related to Z 2 topological 
insulators [SVI IFMP2| . However, our methods can be adapted in order to deal with the two 
dimensional fermionic case ([GMT], in preparation). 

The discussion above is about the case when no magnetic fields are present. While the 
existence of Wannier bases for real and periodic Schrodinger operators can be reduced to the study 
of the existence of smooth and periodic Bloch bases, for non-periodic systems the situation is more 
complicated. Nevertheless, for a large class of non-periodic perturbations of periodic systems it 
can be proved that localized bases still exist, see |NN2) : the main idea behind it is based on a 
‘continuity’ argument which interpolates between the periodic and non-periodic operators. 

Our second main result gives the construction of bases of exponentially localized functions for 
the case when a weak, globally bounded magnetic field is added, taking for granted the existence 
of such bases for the zero magnetic field case. 

The content of the paper is as follows. In par aeranh l 1.2 1 we outline how one arrives at the family 
P(k) starting from a periodic tight-binding Hamiltonian/ Schrddinger operator. In paragraph 11.31 
we formulate the main results. Sections [2] and [3] are devoted to proofs. 
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1.2 Prom time reversal symmetric Hamiltonians to periodic P(k) satis¬ 
fying conjugation symmetry 

Let S C be a finite set containing D points and let 

A := S + Z'* 

be the discrete configuration space. We assume that every point y ^ A can be uniquely written 
as: 

y = y + 7, yeS, 7 g Z"*. 

Consider the Hilbert space and let Hq be a bounded self-adjoint operator acting on this 

space, which commutes with the translations acting on Z"^. It is uniquely determined by the 
’matrix elements’ 


-ffo(y,7;y',7') = iio(y,7-7';y',0), Vy,y'GS, V7,7 'gZ''. 

Consider the unitary Bloch-Floquet-Gelfand transform 

;2(a) {G^){y; k) := ^ + 7 ) G A"([-1/2, 1/2]^; C^). 

7eZ‘^ 


We have 


GHoG*= h(k)dk, h(k)GC{C^), h(j/,y';k) := y e2-k-^ilo(y,7;y',0). 

J[-II2M2V 

The D X D matrix /i(k) will have D eigenvalues (Bloch energies), whose range determine the 
spectrum of Hq. We see that h(k) can be periodically extended to the whole of and that 
the decay of i^o(^/, 7 ;y^ 0 ) as function of 7 tells us how regular /i(k) is. For example, exponential 
localization of Hq gives real analyticity for h. Moreover, if Hq is a real kernel, then we immediately 
obtain: 

= ^(y,y';k) = /i(y,y';-k), h(k)/ = /i(-k) 7 , 

a property which from now on will be called conjugation symmetry. We will see in Lemma 12.51 
that for a self-adjoint matrix, conjugation symmetry is equivalent to: 


*/i(k) = h{-k), k G 


where *h(k) means transposition. 

Assume that the union of the ranges ot N < D eigenvalues Aj(k) of h{-) is separated from the 
range of the other eigenvalues, i.e. there exists a gap in the spectrum of Hq. We can define F’(k) 
to be the Riesz projection associated to these N eigenvalues. P(-) is periodic and has the same 
type of smoothness as h. Moreover, because 

^{{hik)-z)-^} = {h{-k)-z)-\ 

it follows that P(-) has the same conjugation property. 

Consider an orthonormal system of N vectors Sj(k) G which form a basis for Ran(P(k)), 
continuous in k G [—1/2,1/2]^. Then we can define: 

= Wj{y + -f) ■■= [ e"^’"**"'>'Sj(y;k)dk, 1</<V y&B, 7 G Z"*. (1.2) 

J[-l/ 2 ,l/ 2 ]‘^ 

It is known that the set obtained by translating them as follows: 

{wj{- - 7 ) : 7 G 1 < J < iV} C l^iA) 
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forms an orthonormal basis for the range of the spectral projection of Hq corresponding to the 
part of the spectrum given by 

N 

o-Q := Ran Aj(-). 

A similar problem can be formulated for continuous operators. Let R be a Z'^-periodic, real 
and bounded potential. The operator H = —A + V has purely absolutely continuous spectrum. 
Let O := [—1/2,1/2]'^. Consider the Bloch transform 

^ (F«')(y;k) := ^ ^'(y+ 7) £ L^{n-L^{n)). 

7eZ‘i 


We have (in the formula below, k-BC is a short-cut for k-dependent boundary conditions f{y+'y) = 
when both y and y -I- 7 belong to the boundary of fl): 

FHF* = f /ikdk, /ik = — A-I-R with k—BC in 

Jn 

Each fiber Hamiltonian /ik has purely discrete spectrum. Assume that the range of N eigenvalues 
form an isolated spectral island (Tq- Using Combes-Thomas exponential estimates and elliptic 
regularity (see Proposition 3.1 in [CN2] L one can show that the Riesz projection H corresponding 
to tTo has a real integral kernel n(x, x') which is jointly continuous and exponentially localized 
around the diagonal, i.e. there exists a > 0 such that |n(x, x')| < L Then 


FUF* 



^k(y,y') = E e^^^^ m{y + r,]/) mL\n). 

7eZ'i 


Again, we see that Pk has the conjugation symmetry, is real analytic and Z'^-periodic in k. Note 
that (11.21) also makes sense in the continuous case, with the only difference that y € H. 

In the continuous case one can also consider the Bloch-Zak transform 


p 2 (R'i) 9 (P 2 ^')(y;k) := ^ «'(y-f 7 ) S {VL-{n)). 


At fixed k and for G C'(j”(]R‘^), the function (Pz'I')(-; k) is C^iVl) and periodic. We have: 


/*© 

FzHF^ = / h^dk, (ik = (—fV -I- 27rk)^ -|- V with periodic — BC in L^(fl). 

Jq 

The Riesz projection H is decomposed as: 

FzHF* = Tnkdk, nk(y,y') = ^ e^”'^<y+^-^My + r,]/) inL^in). 

7eZ‘i 

We see that the range of Hk consists of H-periodic functions in y, but Hk is no longer periodic in 
k. Instead, denoting by t\ the unitary operator acting on L^(H) given by {T\f){y) = 
we have the identity: 

Hk+A = rAnkr;^\ VAeZ''. 


If the Sj’s are eigenvectors of /i(k), then the corresponding wj’s (see (11.21) ') are called Wannier 
functions. In general, the Wannier functions are just square integrable, but in many physical 
applications it is important to construct better spatially localized Wannier functions. For example, 
if we can simultaneously choose the Sj’s to be periodic and real-analytic, then a standard Paley- 
Wiener argument shows that the Wannier functions are exponentially localized. Moreover, if we 
also have Sj(k) = Sj(—k), then the Wannier functions are real. 

If the Sj’s are not eigenvectors for the fibre Hamiltonian but they form an orthonormal basis 
for Ran(P(k)), we can still define the corresponding Wj’s. They are called composite Wannier 
functions. 
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1.3 Construction of composite Wannier functions 

If 9 k !->• P(k) G C{'H) is a projection valued continuous map in some separable complex 
Hilbert space 71, then the dimension of Ran(P(k)) is constant. We list four basic problems. 

Problem 1.1. Let d>\ and let 9 k P(k) G be an orthogonal projeetion valued map 

which is -periodic and belongs to C'”(R‘^), n > 1. Let dimRan(P(k)) = N < oo. Can one find 
a system of N vectors {Sj(k)}j^]^ which form an orthonormal basis in Ran(P(k)) and which are 
both in C'"(R'^) and if -periodic? 


Problem 1.2. If d > 1, let R"^ 9 k h-> P(k) G P(P) be an orthogonal projection valued map, 
-periodic and jointly analytic in a tubular complex neighborhood I ofMf. Let dimRan(P(k)) = 
N < oo. Can one find a system of N vectors {Sj(k)}j^j^ which are jointly analytic in a (possibly 
smaller than I) tubular complex neighborhood I' ofW^, are if-periodic and form an orthonormal 
basis o/Ran(P(k)) when k G R^? 

Most of the results of this paper will only be valid if some conjugation symmetries hold. Here 
is a list of definitions: 

Definition 1.3. We say that a family Hf) of linear operators acting on a complex Hilbert space 
H has the CS property if there exists an anti-unitary operator K such that 

/C = Id, KH(k)K = P(-k), k G R'^. (1.3) 

The most common example of such a RT in P{1) is the complex conjugation Kf = f: 

PMf = P(-k)7, kGR'". (1.4) 

In fact, up to a change of basis, every such K is given by a complex conjugation. 

Definition 1.4. We say that a family of invertible operators U{-) aeting on a eomplex Hilbert 
space H has the CS' property if there exists an anti-unitary operator K such that 

= Id, KU{\f)K = [P(-k)]-7 k G R”^. (1.5) 

The third problem is the following: 

Problem 1.5. Assume that P(-) has the CS property. Can we construct a basis as in the previous 
two problems, which also obeys the property A"E!j(k) = E!j(—k)? 


Finally, let us formulate a more general problem as in [Pi IFMP] . motivated by the Bloch-Zak 
transform from the continuous case. Assume that we have some separable complex Hilbert space 
% with a conjugation K as in (11.311 . and a family of unitary operators t\ with X (^lf such that 

To = Id, = ta+^(= t^ta), Ta = Tf^{= t_a), Ktx = T-xK. (1.6) 

Problem 1.6. If d > 1, let R"^ 9 k n(k) be an orthogonal projection valued map in H. such 
that: 

n(k + A) = TAn(k)T;, Kii(k)K = n(-k), vk g va g i‘^. 

Assume that the map is jointly analytic in a tubular complex neighborhood I ofRf. Assume that 
the dimension of Ran(n(k)) equals N < oo. Can one find a system of N vectors {'I'j(k)}j^j^ 
which are jointly analytic in a (possibly smaller than I) tubular eomplex neighborhood I' of 
which form an orthonormal basis o/Ran(n(k)) when k G R'^, obey rA4'j(k) = 'kj(k +A) and 
Rr4'j(k) = 4'j(—k) for all k G R'^ and X & if ? 
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Remark 1.7. We will show in Subsection (EH) that this more general problem can be reduced 


to the periodic case. However, other more subtle properties like the existence of a canonical Berry 
connection only appear in the Zak picture, see [FCGj for details. 

Here is our first main result. 

Theorem 1.8. The following statements hold true: 

(i) Problem can be reduced to the first three; see Subsection \2.1i 

(ii) All problems can be reduced to finite dimensional Hilbert spaces; see Subsection \2.^ 

(hi) Assume that the family P(k) is real analytic. If we can construct a continuous solution to 
our problems, then it can be made real analytic, preserving all its other properties; see Lemmas 

[QandrgH- 

(iv) If d =1, a solution exists even in the absence of the CS property; see (12.141) and (I2.17I) .’ 

(v) If d = 2, a solution can be constructed if the CS property holds true. The main technical result 
is Provosition \2. 1 ^ 

Remark 1.9. Let us try to explain very roughly the main ideas behind the proof of (v). The 
construction is inductive in d; using (iv) we can construct some basis vectors 'I’j{ki,k 2 ) which 
are continuous in k but periodic only in ^ 2 - Due to the periodicity of the projection P(k), the 
basis at ki = 1/2 is related to the basis at fci = —1/2 by a matching N x N unitary matrix 
/3(fc2) which is continuous and periodic in ^ 2 , and whose matrix elements have the symmetry 
[l3{k2)]mn = [Id{—k 2 )]nm- The main question is how to rotate each 'i>j{ki,k 2 ) using I3{k2) such 
that the new vectors to be equal at fei = ±1/2. This is done in Proposition l2.16l a result interesting 
in itself. Also, a related question is finding sufficient conditions for a matrix like /?(fc 2 ) to admit a 
periodic and continuous logarithm in k 2 . The main obstacle to the existence of such a logarithm 
is the possible crossing of eigenvalues of /3(fc2) when ^2 varies, see ()2.25|1 for a generic example. 
One of the crucial ingredients in order to circumvent the crossing problem is the Analytic Rellich 
Theorem [R-S4] which allows us to continuously follow the eigenvalues through crossings, also 
preserving periodicity in k 2 . 

Remark 1.10. The lack of a Rellich Theorem in more than one dimension is also the main 
reason for which our proof cannot be generalized to d = 3. We conjecture though that a variant of 
Proposition 12.161 remains true if d = 3, but the proof becomes much more technical and seems to 
be necessary to appeal to ‘avoided crossing’ techniques involving Sard’s lemma. This construction 
will be done elsewhere. 

1.4 The non-zero magnetic field case 

Our second main result of this paper concerns the question as to whether the existence of Wannier- 
type localized bases is stable under the application of a (weak) non-decaying magnetic field. For 
simplicity, we only work in The setting is as follows. We assume that there exists a background 
magnetic field, orthogonal to the plane, whose only non-zero component is denoted with Bq{x.) 
and obeys: 


||d?o||ci(R2) < OO. 

We can always construct its associate transverse gauge: 


(1.7) 



( 1 . 8 ) 


Let F be a bounded, real scalar potential and consider the Hamiltonian 


Ho = (-^V - Aof + F. 


(1.9) 
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Assume that the Hamiltonian Hq has an isolated spectral island o-q and let Pq be its corresponding 
spectral projection. We also assume that the subspace Ran(Po) has an (exponentially localized) 
generalized Wannier basis; let us explain in detail the meaning of this. Let F C be an infinite 
set of points such that 

inf{|| 7 - 7 '|| : 7 , 7 ' e F, 7 ^ 7 '} > 0 . 

The above condition implies that F is countable; in particular, F can be a periodic lattice. Let 
N < 00 and assume that the set of functions 


G Ran(Po), 1 < J < IV, 7 G F 

form an orthonormal basis of Ran(Po). We also assume that there exist a > 0 and M < 00 such 
that 


sup [ |wj-^(x)|2e2“ll^”'^lldx < M. 

i,7 


Now let us consider a magnetic field perturbation given by P(x) such that 

ll^llcMiK^) < 1, A(x) := / sB{s:>^)ds (x 2 ,-:ci), 

^0 

and introduce the perturbed Hamiltonian: 

Hb := (-iV - Ao - bAf + F, 6 G R. 


( 1 . 10 ) 


( 1 . 11 ) 


( 1 . 12 ) 


Below we list a few known facts about the resolvent of Jib, see for details |N3l O ICNj : 


• the resolvent (Hb — z)~^ has a Schwartz kernel denoted by (lib — z)~^(x,x') which is jointly 
continuous outside the diagonal, has a singularity of the type — ln(||x — x'||) near the di¬ 
agonal, and there exist some C(z) < 00 and a(z) > 0 such that |(iFb — z)“^(x,x')| < 
C(z)e-“(^)ll^-^'ll if ||x - x'll > 1 . 




Fix a compact K C p(-ffo)- Then there exist bo > 0, a > 0 
0 < 1^1 < bo we have K C p{Hb) and uniformly in x 7 ^ xh 


sup 

z^K 


(Hb - z)-\x,x') - - z)-7x,x') 


and C < 00 such that for every 
< C | 6 | (1.13) 


where 


</>( 


X, x') := f A{x'+ s{x — x')) ■ (x — x')ds =—4>(x',x). 

Jo 


(1.14) 


In particular, if | 6 | is small enough then Hb has an isolated spectral island ab close to CTq, see also 
[Nl] for a different proof of this fact. Denote by Pb the corresponding projection. Here is the 
second main result of our paper: 


Theorem 1.11. 

(i). If \b\ is sufficiently small, then there exist a constant C < oo, a family of continuous functions 
o,nd a > 0 satisfying 

(1.15) 

such that {Sy 7 ,f,}je{i,...,Ar}, 7 gr o-n orthonormal basis o/Ran(Pb) and 

sup < (j 1 ^ 1 ^ (1-16) 

7J' 


(ii). Moreover, ifV is a periodic lattice, Bq = 0, V is T-periodic, B is constant and 

Wj^fx) = Wjp{x - 7 ) = Wj,fx), (1.17) 
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then 


(1.18) 


and 






(1.19) 


Remark 1.12. When = 1 (the case of a simple band), Theorem 11.111 111 was proved in [N2] . 
In the current manuscript we use a strategy which is related to the one of [N2] ; what we add new 
is a significant generalization of the results and a substantial simplification of the proof due to the 
regularized magnetic perturbation theory as developed in [Cl IN3] . In fact, the estimate (11.131) is 
the main new ingredient of the proof of Theorem 11.111 

Theorem ll.lll is interesting in itself but it also plays an important role in the so-called Peierls- 
Onsager substitution at small magnetic fields [Pel ILul IPSTl IFT] . These aspects will be treated 
elsewhere. 


2 Proof of Theorem 11.8 

2.1 Reduction of Problem 11.61 to the first three 

We first show that there exists a family of unitary operators itk with k € such that u\ = tx 
if A G the map 9 k Uk G -C('H) is entire, and: 

uo = Id, UkUk' = Uk+k'(= Uk'Uk), Uk = 'Uk^(= w-k), Kui^ = W-klf = u^^K, Vk, k' G 

( 2 . 1 ) 

Denote by fj the vectors of the standard basis in The operator is unitary hence the 
spectral theorem implies that it can be written as 


rf, = f e^^dEM)- 

J ( —7r,7r] 

The operator Mj = (j) dEj{(j)) is bounded, self-adjoint, and r/^ = 

If F is any 27r-periodic smooth function we define F(e**) = F{t) and we have 


Firf,) = Y.F{m){Tf,r, /'(m) 


1 

27r 





( 2 . 2 ) 


(2.3) 


We have KtJ^K = (T/j)”™. If F is real, then F{m) = F{—m). This leads to KF{Tf^)K = F(Tf^) 
for all 27r-periodic smooth and real functions F. By a limiting argument we conclude that the 
same remains true for the spectral measure of hence KMjK = Mj. Moreover, because the 
set {Tfj consists of commuting operators, the same is true for their spectral measures and for 
Now if k = [ki,...,kd] G R'^ we define 

One can verify that (EH is obeyed and because the generators Mj are bounded, the map is also 
entire. Now if we define P(k) := uj^^n(k)itk we see that P{-) is Z'^-periodic, has the same smooth¬ 
ness properties as !!(•) and KP{\i)K = P(—k). Assuming that we have a solution {E!j(k)}^;^ for 
Problem ll.51 then we can easily check that l'j(k) := UkSj(k) is a solution for Problem 11.61 


2.2 Reduction to finite dimensional Hilbert spaces 

Lemma 2.1. Fix d and assume that Prohlem M. 1\ can be solved if the complex Hilbert spaee has a 
finite dimension. Then it can also be solved in any infinitely dimensional separable complex Hilbert 
spaee. 
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Proof. Let L > 1 be an integer and consider a net of (2L + 1)"^ points 

ATl :={k„ G [-l/2,l/2]'': k„ := [ni/(2L),... nd/(2L)], -L<nj<L, G Z, l<j<d}. 

Since P{-) is continuous and Z'^ periodic, then for every e > 0 there exists large enough 
such that for every k G there exists some k„ G A/l^ such that ||-P(k) —P(kn)|| < e. Now define 
the finite dimensional subspace (to be understood as a sum of linear subspaces) 

iJe := ^ Ran(P(k„)), 

and denote by lie its corresponding orthogonal projection. Clearly, neP(k„) = P(k„)ne = P(kn) 
for all k„. Define Pe(k) := neP(k)ne. For every k G we have (the choice of k„ below is such 
that ||P(k) - P(k„)|| < e): 

Pe(k) - P(k) = ne{P(k) - P(k„)}ne + P(k„) - P(k) 

hence ||Pe(k) — P(k)|| < 2e. It means that Pe(k) is close to a projection. If e is small enough, 
then the operator: 


Pe{k) :=-^ f (Pe(k) - z) ^dz 

i-Kl 7|z-1| = 1/2 

is an orthogonal projection whose range is contained in (to see this use (Pe(k) — z)~^ + z~^ = 
2 :“^Pe(k)(Pe(k) — z)~^ in the above formula). Using the resolvent formula 

P,(k) - P(k) = ^ [ (P,(k) - z)-\P,{k) - P(k))(P(k) - z)-^dz 

we obtain ||P(;(k) — P(k)|| < const e. The dimension of P(k) is constant and equal to N, thus the 
dimension of Pe(k) also equals N for all k, provided e is small enough. This is true because we 
can construct the Sz.-Nagy unitary which intertwines between Pe(k) and P(k): 

P,(k) := {P(k)P,(k) + (Id - P(k))(Id - A(k))}{Id - (P,(k) - P(k))2}-i/2^ (2.4) 

with P(k)P,(k) = P,(k)A(k). 

To summarize: we constructed a projection valued map Pe(-) with dimension N which enjoys 
the same smoothness and periodicity properties as P{-), and lives in the finite dimensional space 
H^. Our assumption was that we can solve Problem 11.11 in finite dimensional Hilbert spaces. 
Hence let ,:(k)}j^]^ be an orthonormal basis of Ran(Pe(k)) consisting of smooth and periodic 
vectors. Then a solution to our problem in the infinite dimensional Hilbert space will be given by 
the vectors: 

vl/,(k) := P,(k)^,-,(k). 

□ 


Lemma [2.11 showed how to reduce the problem to a finite dimensional Hilbert space. We will 
now show that the conjugation symmetry is also preserved by the previous construction. 

Lemma 2.2. We refer to the quantities defined in Lemma \2.1[ Assume that P(-) has the CS 
property. Then both Pe(-) and Pe(') have the same property. Moreover, i/iL\E'j.£(k) = 4'j^(—k), 
the same property holds for (k). 

Proof. First we prove that the subspace Pg is invariant with respect to the action of K, i.e. H^ 
commutes with K. We constructed the grid A/l^ so that both k„ and —k„ belong to it for all 
n. Due to the CS property we have P(—k„)iL = KPfkn), which shows that iLRan(P(k„)) C 
Ran(P(—k„)), hence K{Hf) C He. Using = Id we also obtain He C K{He), thus K{He) = He. 
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Now let f G and g € be arbitrary. Then using the anti-unitarity of K together with 
Kf G we obtain: 

(Kglf) = {Kf\g) = 0. 

This shows that K{H^) C and in fact K{H^) = . Then: 

n,/-f (Id - n,) = (Id - noit'n, = o. 

This implies that li^K = KIl^. It follows that Pe(k) = n£P(k)ne has the CS property. We see 
that 

(P,(k) - z)-'^K'i/ = K{P^{-k) - z)-^^ 

which after integration leads to the conclusion that Pe(k) has the CS symmetry. 

The last thing we need to prove is that the Sz.-Nagy unitary C/e(k) also has CS symmetry. This 
is immediately implied by KP{k)K = P(—k) and KP^{'k)K = Pe(—k), which ends the proof. 

□ 


It is possible to construct an orthonormal basis in such that Kfs = fs- If 'I' = ''«/« 

then K'i> = ^ becomes just a complex conjugation. Thus we can assume without 

loss of generality that the Hilbert space is with 2 < D < oo and K is complex conjugation. 
In this case, all operators are characterized by finite dimensional matrices. 

2.3 Lifting continuity to real analyticity 

The next lemma states that if the map P(-) is real analytic, then a solution for Problem 1 1.1 1 can 
always be turned into a solution for Problem ll.2l 

Lemma 2.3. Let 9 k n> P(k) G C{C^) be a -periodic, orthogonal projection valued map 
with dim(RanP(k)) = N < D for all k G K."^. Assume that all the components {Pmn{-)}i<m,n<D 
in the standard basis of have a -periodic, bounded jointly analytic extension to a tubular 
complex neighborhood I of K'^. Assume that {'I'j(k)}j^j^ is an orthonormal basis of RanP(k) 
for all k G Z‘^ -periodic, and continuous. Then there exists holomorphic on a 

tubular complex neighborhood R'^ G I' C I, Z‘^-periodic, which is an orthonormal basis o/RanP(k) 
if kG 

Proof. The function g(k) = + ^j)~^ obeys g(k) = g(—k), it is analytic on the strip 

{zgC^ : \lm{zj)\ <1, j G {!,..., d}}, and /jjdg(k)dk = 1. If (5 > 0 we define an approximation 
of the Dirac distribution ^^(k) = S~‘^g{6~^k). Let us consider 

Vl/,, 5 (k)= [ 55 (k-k')vl',(k')dk', 


which defines a D dimensional vector which is Z'^-periodic. 

Let us show that admits an analytic extension to a strip 

{z = (zi, ...Zd) G : |Im(zj)| < d}. 

Due to periodicity, it is enough to find an extension to a bounded open neighbourhood O of [0, 
of the form 

O = Oi X ... X Od, Oj = {zj = Xj + iyj G C : —2 < Xj < 2, ly^j < S' < d}. 

For every z = {xi -\- iyi, ■■■, Xd + iyd) G O and every k' G R'^ we have the estimate: 

|«(z - k')l < 
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Also: 


gs{k - k') = 


1 f f 9s{C - k') 

{2Tny Jqq^ Jqq^ nLi(Ci - kj) 


dC, k € [0,1]^ k' e 


Because the 'hj’s are bounded on 


the Fubini Theorem and the Cauchy Integral Formula 
allow one to analytically extend each dtj , 5 (-) to O, and by periodicity, to a strip around 
Let us now continue the construction. A standard approximation argument implies: 


lim i sup H^'j,^(k) - «'j(k)|| 1=0. 

Denote by <l)j_ 5 (k) := P(k)'l>j_ 5 (k). Then for a given 0 < e <C 1 there exists S small enough such 
that: 

max sup \\^j,sik) - «'i(k)|| < e. 
keR'i 

Define the selfadjoint N x N Gram-Schmidt matrix 

/ifej(k) = ($j, 5 (k),$fe, 5 (k)) = [ gs{k-y)gs{k-y'){P{k)-^j{y),'iik{y'))dydy'. (2.5) 

jR2d 


Since {'I'j(k)}j^;^ is an orthonormal system, by choosing S small enough we can insure that 
||/i(k) — Id|| < 1/2 uniformly for k € Moreover, h{-) is bounded, Z'^-periodic and real analytic, 
and the same holds true for Then the vectors 


N 

:= 

fc=i 


1/2 


kj 


‘hfc,(5(k) 


form a basis with all the required properties. 


□ 


2.4 Preservation of 

Lemma 2.4. Let T’(k), and as constructed in Lemma \2.S\ Assume that 

P(k) has the CS property, and that dij/k) = 4'j(—k) for all j. Then Sj(k) = E!j(— k) on for 
all j. 

Proof. We will refer to the already defined objects in the proof of Lemma [2.31 First, due to the 
symmetry of g it is easy to check that 'I'j_ 5 (k) = k). This implies that <I)j, 5 (k) = k). 

This leads to the identity hjkik) = hjk{—k). Since /i(k)~^/^ can be written as a power series in 
h{k) — Id, we have that [h{k)~'^/'^]kj = \h{—k)~^/‘^]kj. □ 

2.5 Some facts about matrices having CS or CS' symmetry 

Lemma 2.5. Let 9 denote complex conjugation and consider a self-adjoint matrix family a(k) £ 
C{C^) having CS symmetry, i.e. 9a{k)9 = «(—k). Then: 

[Q;(k)]m„ = [a{-k)]nm, cr(a(k)) = tT(a(-k)). (2.6) 

Proof. We know that [Q:(k)]j„n = [Q:(k)]„j„ from self-adjointness and [a(k)]mr!, = [ci;(— k)]mn from 
CS symmetry. The symmetry of the spectrum follows from the identity 

det(A — 2 Id) = det(*A — 2 ; Id). 


□ 
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Lemma 2.6. Let 9 denote the complex conjugation and consider a unitary matrix family /3(k) £ 
C{C^) with the CS' symmetry, i.e. 0/3(k)0 = /3(—k)“^ for all k. Then: 

[/3(k)]m„ = [/3(-k)]„™, cr(/3(k)) = cr(/3(-k)). (2.7) 

Proof. We know that [/3(-k)-i]m„ = [/3(-k)]„m from unitarity and [/3(-k)-i]m„ = [/3(k)]m„ from 
the CS' symmetry. The proof of the symmetry of the spectrum is the same as in the previous 
lemma. □ 


Remark 2.7. We note the fact that both the CS property for self-adjoint matrices and the CS' 
property for unitary matrices boil down to the same identity, only involving transposition: 

*a(k) = a(—k). 


Lemma 2.8. Let k = [ki, /C 2 ] £ and let Q:(k) be self-adjoint with CS symmetry. Then the uni¬ 
tary family /3(k) := has CS symmetry, while the family 7 (k) := has CS' symmetry. 


Proof. We have 
In a similar way: 


_ ^-iki9a{k)0 _ 

e-f{k)9 = = [7(-k)]-i. 


□ 


2.6 The case d = 1 

Let us start by stating without proof a well-known fundamental lemma essentially due to Kato, 
formulated in a form which will be convenient for us in what follows. 

Lemma 2.9. Let Q{x) = Q*{x) = Q‘^{x) be a C^-norm differentiable family of orthogonal projec¬ 
tions. Define 


lC{x) := i(Id - 2Q{x))Q'{x) = i[Q'{x), Q(x)], lC{x) = lC*[x). 
Given y £'S., one can define a family of unitary operators A(x, y) such that 

idxA{x,y) =lC{x)A{x,y), A{y,y) = ld, 


( 2 . 8 ) 


(2.9) 


given by: 


A{x,y) f dsi f ds 2 -.- f dsnJC{si)..X{sn)- 

~J >1 dy Jy Jy 


Furthermore, we have: 

idyA{x,y) = -A{x,y)H{y), [A[x,y)]~^ = A{y,x), A[xi,X 2 )A{x 2 ,xf) = A{xi,xf). 
The key intertwining identity is: 

Q{x) = A{x,xo)Qixo)A{xQ,x). 


( 2 . 10 ) 

( 2 . 11 ) 

( 2 . 12 ) 


For fixed xq £ K, if either Q{-) £ (^"(R) with n > 1 or Q{-) admits an analytic extension to the 
strip := {z £ C : |Im( 2 ;)| < a}, then the same holds true for A(-,xo) and A{xtj, •). 
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In the case in which the above family is Z-periodic, i.e. Q(x + 1) = Q{x) for all x S M, 
then /C(x) is Z-periodic and we observe that A{x -b l,xo -b 1) must equal A(x,Xo) (see (12.91) 1. 
Consider the unitary operator A(1,0). Reasoning as in (12.21) . the spectral theorem provides us 
with a self-adjoint operator M with spectrum in (—7r,7r] such that bl(l,0) = Also, due to 

the intertwining relation Q(1)A(1,0) = A(l,0)(5(0) and because <5(1) = Q(0) due to periodicity, 
we conclude that A(1,0) commutes with <5(0) and so does M. 

Define the unitary operator: 

U{k) := A{k, 0)6“*'=^, fc e K. (2.13) 

One can verify that !/(•) is Z-periodic: 

U{k -b 1) = A{k + 1,= A{k + 1,1)A(1,= A{k, = C/(/c) 

where we used A{k -b 1,1) = A{k, 0). 

Now let us consider Problem 11.11 with d = 1. Consider any orthonormal basis {^'j(0)}^i in 
the range of P{0). Construct U{k) as above starting from P{k). We will now prove that the 
vectors: 

Sj(fc) := C/(fc)«'j(0), fceK, l<j<N (2.14) 

give a solution for Problem ll.il These vectors are clearly periodic and orthogonal, thus we only 
need to prove that they live in the range of P{k). But since M commutes with P(0) and due to 
the intertwining relation P{k)A{k, 0) = A{k, O)P(O), it follows that P{k)U{k) = U{k)P(0) and we 
are done. 

We note that the above construction has not used a possible CS property of P(-). If such a 
property does hold, i.e. besides smoothness and periodicity we also have 9P{k)9 = P{—k) for all 
k, then let us show that we can construct a basis which solves Problem ll.5l 

Since 9P{0) = P(O)0, we have that 0'I'j(O) belongs to the range of P(0), hence the real 
vectors ^(th^^O) -b04'j(O)) and A(\]/^.(o) — 0vl/^(o)) have the same property. After a Gram-Schmidt 
procedure we obtain a real orthonormal basis of Ran(P(0)), denoted by {Sj(0)}j^]^. Thus the 
only thing we need to check is that U{k) has the CS property, i.e. 9U(k)9 = U{—k). 

First, due to the identity P'{—k) = —9P'{k)9 we have 9K.(k)9 = /C(—fc), thus 9A(k,0)9 = 
A(—fc,0) because they solve the same initial value problem. 

Second, taking fc = 1 in the last identity we get 9A{\,Q)9 = A(—1,0) = A(0,1) = [A(1,0)]“^ 
i.e. 

0C^9 = e-^^, (2.15) 

where M = M* commutes with P(0) and has its spectrum in (—tt, tt]. Since the underlying Hilbert 
space is C^, all operators are represented hy DxD matrices, hence M has discrete eigenvalues and 
can be written as M = J2-Tr<<i><Tr where H^ are spectral orthogonal projections corresponding 
to different eigenvalues. From USTSl) we immediately obtain that = H^ for all eigenvalues, 
hence 9M9 = M. Thus: 

Q^-^kMQ ^ ^^keMe ^ ^ikM ^ fceR, 

hence 9U{k)9 = U{—k) and the basis vectors we are looking for are given by Sj(/c) := P(A:)Sj(0). 

Remark 2.10. Define ’^mik) := A(fc, O)Sm(O) and consider Mjk ■= (MSfc(0)|Sj (0)), i.e. the 
matrix elements of P(0)MP(0) in the real basis S(0). Using the fact that M commutes with 
P(0) = P(l) leads to 

N 

^„(fc + 1) = A(fc, 0)P(0)e*'^P(0)S„(0) = '^[P^]nm'fn{k) (2.16) 


S^(fc) = P(fc)S,„(0) = A{k, 0)P{0)e-^^^E^i0) = ^[e-*^'^]„„vi;„(fc). (2.17) 

n—1 

In other words, we have found a unitary, fe-dependent change of basis in the range of P(fc), which 
transforms the 'I'(fc)’s into a periodic basis also preserving the CS symmetry and the regularity 
of the original projection. 
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2.7 The case d = 2 


From now on (if not otherwise stated) we assume that the underlying Hilbert space is and 

the family of projections P(-) has the CS property, i.e. 9P{h)9 = P(—k) for all k S The 

idea is to proceed recursively. Consider P(0, /c 2 )- Using the result from d = 1 we can construct an 
orthonormal basis {Sj(0, ^ 2 )}^! of RanP(0, ^ 2 ) such that 

Sj(0,fc2 + l) = Sj(0,fc2). 9Ej{0,k2)9 = Ej{0,-k2), Vfe e R. (2.18) 

Now keeping ^2 fixed, consider the orthogonal projection family n(-) := P{-,k 2 ). We can apply 
Lemma |2.9I with Q replaced by H and obtain a family of unitary operators denoted by Ak^ix,!/) 
generated by the self-adjoint operators 


lCk^{ki) := t[5iP(fci,fc2),P(fci,fc2)]. 

Reasoning as in the case d = 1 we obtain the intertwining relation; 


P(/ci,fc2)Afe,(fci,0) =Hfc,(fci,0)P(0,fc2), (2.19) 

and 

d/Cfc,(fci)d = /C_fc,(-fci), dHfc,(A:i,0)d = H_fc,(-fei,0). (2.20) 

We also note the periodicity in k 2 '. 

/Cfc,+i(fci) =/Cfc,(fci), Hfc,+i(fci,0) =Hfc,(fci,0). (2.21) 

Due to the intertwining relation (I2.19|) we have that the set of vectors 

5'm(k) := Hfe2(fci,0)Sm(0,fc2) 

form an orthonormal basis in RanP(k), having all the right properties but one: they are not 
periodic in ki. Using Ak^{ki -|- 1,1) = Ak^{ki,9i) we observe that: 

^'m(fcl + 1 ,^ 2 ) = Hfc2(fci,0)Hfc2(l,0)Sm(0,/c2). 

From the intertwinine relation (12.191) and due to Pfl.fe9l = P(0.kA we have that AkAhO) 
commutes with P(0, /C 2 ) hence: 

N 

'fmikl + 1, k2) = Y.^Ak, (1,0)S„(0, fc2)|S„(0, fc2))^n(fci, ^ 2 )- (2.22) 

n—1 

The N X N matrix family defined by 

Pnm{k2) := {Ak, (1, 0)S„(0, fc2)|S„(0, k2)) (2.23) 

is unitary in and we have: 

N 

^'m(fcl + 1, ^ 2 ) = y] dnm{k2) ^'«(fcl, ^ 2 )- (2.24) 

n—1 

This identity closely resembles (j2.16|) but the important difference is that the matrix /3 depends 
on /c 2 . The important question we have to answer is whether we can find a k-dependent unitary 
rotation, like we did in p.l7L which transforms the 'k’s into a Z^-periodic basis, preserving at 
least its continuity and CS symmetry. We note that d{k 2 ) is similar to the ‘obstruction unitary’ 
Uahs(k 2 ) introduced in [FMPj . 

Lemma 2.11. The N x N matrix family defined hy (12.2311 is unitary, h-periodie, inherits the 
regularity properties o/P(0, fc 2 ) and has the CS' symmetry. 
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Proof. We only prove the statement related to the CS' symmetry. Since 


( 1 , 0 ) = 0Ak, (- 1 , 0)9 = OAk, ( 0 , 1)0 = ( 1 , 0)]*0 

and because 02 ^( 0 , —^ 2 ) = Sm( 0 , ^ 2 ) we have: 

/lnm(-fc2) = (0[Afc,(l,O)]*S^(O,fc2)|0S„(O,fc2)) = (S„ (0, ^ 2 ) | [Afc, (1, 0)] (0, ^ 2 )) = /3™„(fc2), 

which according to Lemma 12.61 is equivalent with the CS' symmetry. □ 


Remark 2.12. Let us assume that there exists an N x N self-adjoint matrix family h{k 2 ) which 
is continuous, Z-periodic, with the CS symmetry, such that f3{k2) = Note that this was 

the case when d = 1 (there h(k 2 ) = M, see (|2.16p L Then the matrix has the CS' 

symmetry (see Lemma [TS]) and the vectors 

N 

s™(k) 

n—1 

give the periodic, continuous and CS symmetric basis we are looking for. Indeed, in this case we 
have: 

N N 

S^(fci + l,fc2) = = S^(k). 

n=l j=l 


Unfortunately, finding a ‘good’ logarithm for / 3 (fc 2 ) (i.e. a matrix h(k 2 ) which is self-adjoint, 
periodic, continuous, CS symmetric and / 3 (/c 2 ) = is not always possible, as it can be seen 

from the following example. Let 


m 


cos(27rA:) sin(27rA:) 
— sin(27rfc) cos(27rfc) 


cos(27rfc)Id -I- i sin(27rfc)cr2, fc G R 


(2.25) 


be a family of unitary matrices, real analytic, Z-periodic, and which obeys the CS' symmetry. 
Note that I3{k) has degenerate eigenvalues at fc = 0 and k = ±1/2, and the range of its spectrum 
covers the whole unit circle. 

Let us show that one cannot find a continuous self-adjoint family h{k) which is also Z-periodic, 
has the CS symmetry and fi{k) = for all k. Assume the contrary. From det(/l(fc)) = 1 = 

giTr(/i(fc)) have that k i—>■ ^Tr(/i(A:)) is a continuous, integer valued function. Hence there exists 

an integer m such that Tr(h(fc)) = 27rm for all k. Let Fj{k) := then we have 


h{k) 


— ^^ Id + Fj(k)uj = mTrld -|- F(fc) • cr. 
i=i 


Thus 

m = (-l)-e®^W- = (-1)™ (cos(|F(fc)|)Id + • a) . 

Comparing with the definition of /3(fc), both Fi and ±3 must be identically zero. Hence F(fc) = 
[ 0 ,F 2 (fc), 0 ] and |F(fc)| = |F 2 (fc)|- Therefore 


/3(/c) = (—1)™ (cos(F2(fc))ld ± isin(F2(fc))o’2) = cos(F2(fc) ± wTrjld-I-i sin(F2(fc) ± TO7r)o-2- 


We conclude that there exists an integer-valued function n{k) such that F 2 {k)+m'K = 27r/c-|-27rn(fc). 
Since ±2 must be continuous, it follows that n(k) is constant and equals some n. Hence F 2 {k) = 
2'Kk + 7 r( 2 n — m) for all k. In order to make sure that j3pq(k) = f3qp{—k) we must have that 
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F 2 {k) = —F 2 {—k), i.e. m = 2n, (—1)™ = 1 and F 2 {k) = 27r/c. This shows that F 2 cannot be 
periodic. 

Fortunately, the lack of a ‘good’ logarithm can be circumvented by a two step rotation pro¬ 
cedure (see (12.341) 1. where for each step a ‘good’ logarithm can be constructed. In what follows, 
we will discuss four particular situations in which a good logarithm can be constructed and then 
show how they can be combined in order to deal with the general case. 

2.7.1 When N = 1 

Here the unitary matrix P{k 2 ) is just a complex number on the unit circle, which due to the CS' 
property (see (12.71) 1 is also even: /3(fc2) = /3(—^ 2 ). 

Lemma 2.13. There exists a unique real and even function 4'p{k2) = 4'p{~k2), with the same 
smoothness properties as fi, Ij-periodic, (j)p{0) = 0, such that P{k 2 ) = 

Proof. Without loss of generality we may assume that /3(0) = 1. We start by first proving the 
existence of a continuous phase <f>{t) such that /3(t) = e*"^*^*^ and (j){0) = 0. 

The map /3 is uniformly continuous on R, thus there exists (5 > 0 such that 

\P(t) — /3(s)| = |/3(<)//3(s) — 1| < 1/2 whenever |t — s| < S. 

If Ln denotes the principal branch of the natural logarithm, then in a neighbourhood [—<5, <5] we 
can define (j){t) = —zLn(l -|- (/3(t) — 1)) = </>(—t) where /3{t) = ^(0) = 0, ^ is real and 

continuous. 

We can repeat this construction in the interval {6, 2S] by defining 

(j){t) = cj){S) - iLnjl -b {P{t)/l3{S) - 1)]. 

A similar formula holds on [—2^, —5), the extension being continuous at and symmetric. After 
a finite number of steps we can cover any bounded interval of R. Thus we obtain a continuous, 
symmetric, real function (f with ^(0) = 0, not necessarily periodic, such that /3(t) = We 

only need to prove its Z-periodicity. 

We know that /3 is Z-periodic, thus the function 

•= - 2 /^- 

must be both continuous and integer valued. Thus F(t) equals some constant n € Z. But we 
have n = F(—1/2) = 0 because (j){\/2) = </(—1/2). This shows that (f>{t -|- 1) = (j)(t) for all t. The 
regularity properties of (j) are the same as those of /3, which can be seen from the formula 

(j,{t) = (j){to) - iLn[l -b {I3{f)/j3{to) - 1)] 

valid for all t with |t — to| small enough. In particular, if (3 is analytic in a strip containing R, the 
same holds true for (f. 

The uniqueness of such a phase is a consequence of the following fact: if is continuous, 
/>(0) = 0, and = 1 for all t, then (j) must be a constant equal to zero. □ 

2.7.2 When the eigenvalues of /3 never cross 

Here we assume that > 1 and /3(fc) has nondegenerate spectrum for all k, besides being 
continuous, Z-periodic and with CS' symmetry. 

Lemma 2.14. There exists a self-adjoint N x N matrix family h(k) which is continuous and 
Z-periodic, has the CS property and /3{k) = 
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Proof. P is uniformly continuous on [—1/2,1/2]. Given any e > 0 there exists S > 0 such that 
||/3(/c) — P{k')\\ < e whenever \k — k'\ < 6. Because the eigenvalues never cross, there must exist 
some positive distance eg > 0 between any two eigenvalues at a given k. Let Jq > 0 be such that 
\\P{k) — /3(fc')|| < eo/10 whenever \k — k'\ < Jq. 

Because /3(0) has N nondegenerate eigenvalues, we can order them according to their argu¬ 
ments. Thus we have —tt < pi < (j )2 < ■■■ < 4>n < and Aj(0) := We know that 

||/3(fc) — /3(0)|| < eo/10 whenever jfcj < Jq- Regular perturbation theory implies that each set 
a{P{k)) n consists of exactly one point when jfcj < Jq. The operators 

Pjik) ■= tt f {P{k) - z)~'^dz, |fc|<5o, 

J\z-\jm=eo/5 

are the orthogonal spectral projections of P on \k\ < i5o. They are continuous. Moreover, 

Xj{k):=Tr{p{k)P,{k)), |fc|<<5o 

are N continuous eigenvalues. Also, because the spectrum of P is symmetric (see (12.Tl) !. we must 
have Xj{k) = Xj{—k) on jfcj < Sq- Now we can repeat the procedure on [—2(5o, —(^o] U [(5o,2i5o] 
starting from k = ±i5o. 

After a finite number of steps we obtain some continuous eigenvalues Xj{k) = Xj{—k) on 
— l/2<k< 1/2, together with their one dimensional Riesz projections. The disk \z — Xj{k)\ < eo/5 
does not contain any other eigenvalues of P{k) besides Xj{k), uniformly in \k\ < 1/2 and j, and 

P,{k) = ^j {P{k)-z)-^dz, \k\<l/2. 

J|z-A,-(fc)|=eo/5 

Using the fact that *'[{P{k) — z)~^] = {P{—k) — z)~^ (from the CS' property) together with 
Xj{k) = Xj{—k), we have ^Pj{k) = *Pj{—k). Moreover, Pj{—1/2) = Pj{l/2) due to the periodicity 
of p. 

Now any fc £ M can be uniquely written as k = n + k with n £ Z and —1/2 < fc < 1/2. Define 

A,(fc):=A,(,fc), P,{k) := P,{k)- 

Since the spectrum of P is periodic (as a set), the Aj’s are the periodic and continuous eigenvalues 
of P on R, which for n > 0 obey: 


Kj{n -|-fc) = Aj(fc) = Xj{—k) = Aj{—n — fc), 

i.e. they are also even. From Lemma 12.131 it follows that we can find pj periodic, continuous 
and even such that Aj(fc) = _ The corresponding projections, Pj(k), (as given by the Riesz 

formula) are continuous, periodic, and obey 

Pj{n + k) = Pj{k/) = *Pj{—k) = *Pj{—n — fc), Vn > 0. 

Thus h{k) := 4'j{k)Pj{k) obeys all the necessary conditions. 

□ 


2.7.3 When a global Cayley transform exists 

Again we assume TV > 1, P{-) continuous, Z-periodic, and with CS' symmetry. We also assume 
that there exists some po £ (—7r,7r] such that e'“^° never belongs to the spectrum of P{k). Then 
the statement of Lemma [2.141 remains true, although the proof is rather different, as we will show 
in what follows. 

The matrix 7 (fc) := e®^''®“'^“^/3(fc) has the property that —1 never lies in its spectrum. The 
matrix 

s(fc) := T{ld — 7 (fc)}{ld + 7 (fc)}“^ (2.26) 
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is self-adjoint, continuous, periodic and 9s(k)9 = s(—fc), hence it has the CS symmetry. By 
elementary algebra we get: 


7 (fc) = {Id -I- zs(/c)}{Id — *s(fc)} 

Consider f{z) = (1 -I- iz)/{l — iz) defined on a (narrow) strip S containing the (real) spectrum of 
s{k) for all k. Then: 

h{k) := -TT + (j}o — iLn'y{k) = -tt + 4>o + — [ {Lnf{z)){s{k) - z)~ 

JdS 

obeys all the properties. 




dz 


(2.27) 


2.7.4 When /3( ) is real analytic and /3(0) and /3(l/2) have nondegenerate spectrum 


This is the last particular situation in which we prove the existence of a ‘nice’ logarithm with 
the properties listed in Lemma 12.141 The real analyticity of /3(-) allows us to relax the strict 
non-crossing condition we had in Lemma 12.141 now we only demand nondegenerate spectrum at 
k = 0 and k = 1/2. Note that the periodicity of /?(•) implies, in particular, that the spectrum at 
k = —1/2 is also nondegenerate. 

The first important remark is that the spectrum of /?(•) consists of real analytic eigenvalues. Let 
us construct them. Assume that e‘'^° is not in the spectrum of /3(0) and define 7 (fc) := e*^‘"’“'^‘’^/3(fc). 
Then Id -I- 7 (fc) is invertible for k near 0 and we can proceed as in (12.261) and (12.271) and construct 
a self-adjoint matrix h{-) which is analytic in a small open disk containing k = 0 such that 
P{k) = gd-'^+'i’o+Hk)) Q_ yyg ggg label the eigenvalues of /3(fc) as {Aj(/c) = 
choosing for example that at fc = 0 we have — tt < />i(0) < ... < (j)N{0) < tt. 

The second step is to extend these projections and eigenvalues to |A:| < 1/2 so that they 
remain real-analytic. If the eigenvalues of /3 never cross, then the projections and eigenvalues 
we constructed in Lemma 12.141 are the objects we are looking for. But here the eigenvalues can 
cross, and we need to be able to analytically continue the one dimensional projections through 
crossing points. By applying a finite number of times the Analytic Rellich Theorem |R-S4) to 
the ‘local self-adjoint logarithms’ (j2.27|) we can follow analytically each projection Pj{k) and each 
corresponding eigenvalue Xj{-) to |fc| < 1/2. Because (3 is periodic and /3(l/2) = /?(—1/2) are 
non-degenerate, this construction can be extended to the whole of K. 

Near k = 0 the spectrum of (3 is nondegenerate and we can write each Pj(k) as a Riesz integral 
on a /c-independent contour: 


Pj{k) 


i 

27r 


'\z-Xj{0)\^ro 


im 


— z) ^dz, 


hence ^[Pj (fc)] = Pj{—k) near zero, identity which remains true on M through analytic continuation. 

Also, because Aj(fc) = Xj{—k) near fc = 0 (due to the CS' symmetry and the non-degeneracy 
of the spectrum of /3(0)), this equality remains true on K by analytic continuation. 

The only thing we still need to prove is the periodicity of eigenvalues and their eigenprojections. 
We know that Aj(—1/2) = Aj(l/2) = aj from the symmetry, all of them being nondegenerate since 
by assumption /3(l/2)[= /3(—1/2)] has nondegenerate spectrum. Let ro > 0 small enough such 
that \z — cijl < ro does not contain other eigenvalues. Then if fc is small enough we have: 


^dz=-^ f {l3{l/2 + k) — z) ^dz 

J \z—aj\=ra 

due to the periodicity of /3, thus Pj{—1/2 k) = Pj{l/2 + fc) near fc = 0, hence they are equal 
everywhere through real analyticity. This proves that each Pj is periodic. Finally, for small enough 
fc we have: 


J\z-aj\=ro 


(/3(-l/2 + fc)-z) 


Aj (-l/2 + k) = Tr{Pj(-l/2 -t fc)/3(-l/2 -k fc)} = Tr{Pj (l/2 -k fc);9(l/2 -k fc)} = Aj (l/2 -k fc) 

which must hold everywhere so they are periodic. Now we can apply Lemma 12.131 for each Xj 
obtaining Xj{k) = where (pj is continuous, Z-periodic and even. Then the solution is 

h{k) := 
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2.7.5 The general case 


Let us go back to (12.241) . the starting point of our discussion. Remember that the matrix family 
/3(-) has the CS' symmetry and we want to unitarily rotate the so that the new basis stays 

CS symmetric, remains at least continuous and is also periodic. 

Let us rewrite (I2.24p by setting fci equal to —1/2: 

N 

^^(1/2, fcs) = ^[/3(/C2)]r^m ^'n(-l/2, ^ 2 )- (2.28) 

n—1 

We will see in the next lemma that the only relevant interval is fei G [—1/2,1/2]. 

Lemma 2.15. Assume that the vectors {S„(x,/c 2 )}))Li form an orthonormal basis for P(x, ^ 2 ) 
and are continuous on [—1/2,1/2] x K., have the CS symmetry S„(x, fe 2 ) = S„(—x, —/C 2 ) ond 
moreover, S„(—l/2,fc2) = S„(l/2,fc2). Given ki € K., we can uniquely write it as ki = n + x with 
n S Z and x S (—1/2,1/2]. Then the vectors defined by S„(fci, ^2) := S„(a;, k2) are continuous on 
CS symmetric and -periodic. 

Proof. The periodicity is a consequence of the identity S„(x + 71 ,^ 2 ) = S„(a;, ^ 2 ) which is valid 
for all n G Z and x G (—1/2,1/2]. The only thing we still need to prove is the CS symmetry. We 

have: _ 

E„{x + n, k 2 ) = S„(a:, ^ 2 ) = S„(-a;, -^ 2 ) = S„(-a: - n, -/C 2 ). 

□ 


The main technical result needed for the proof of Theorem 1 1.81 is the following: 

Proposition 2.16. There exists a family of unitary matrices u{x, ^ 2 ) continuous on [—1/2,1/2] x 
R and Z-periodic in /c 2 , such that: 

[u(-l/2, fc2)]-^/3(fc2)u(l/2, k2) = Id, (2.29) 

which also has CS symmetry: 

6u{x,k2)0 = u{—x,—k2). (2.30) 

Moreover, the vectors: 


N 

^2) .— ^ ^ [^(^? , ^2) (^’^l) 

m—1 

obey the conditions needed in Lemma \2.15\ 

2.7.6 Proof of Proposition [2.161 

Before anything else, let us relate the proposition to the cases in which we can write /3(fc2) = 
with some ‘nice’ h{k 2 ). In all those cases we can put u{x, ^ 2 ) = and we immediately see 

that (12.291) means no more than: 

^ jd. 

Since a ‘good’ logarithm might not always exist, the new important idea is to realize that one 
can straighten up the initial basis by performing two successive unitary rotations. First, given 
f as in the theorem, we will show that we can always find some ‘good’ (i.e. Z-periodic, with CS 
symmetry and continuous) h{k 2 ) such that 

sup _ id|| < 1/2. (2.32) 

^2 
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Then rotating the ^'-basis with e “^(^= 2 ) we get 'hm(a;, ^ 2 ) := J2n=il^ ^ 2 ) which 

obey: 

N 

$^(1/2, fca) = ^n(-l/2, ^ 2 ). (2.33) 

n—1 

But now the new matching unitary matrix 

/3(fc2) := = Id + [/3(/c2) - 

has the CS' property, is periodic, continuous and always has —1 in its resolvent set. Reasoning as 
in Subsection 12.7.31 (see formula (12.271) 1 we obtain another ‘good’ ( 1 (^ 2 ) such that P{k 2 ) = 

By rotating T'ra(a:, ^ 2 ) with we obtain our E!„(a;, ^ 2 )- This also gives an expression for 

the unitary u: 

u{x, * 2 ) = = Id. (2.34) 

Now let us give the precise statements. 

Lemma 2.17. Assume that we can find M > 1 self-adjoint N x N matrix families Hj{k 2 ) which 
are 1-periodic, continuous and with CS symmetry so that: 

g-tJfi(fc 2 ). ^-iHulk-,) . ^-:LH,{k2) ^ (2 35 ) 

Then u{x,k 2 ) '■= is unitary, h-periodic in the second variable and obeys 

(I2.29|) and (j2.30D . Moreover, the vectors defined in (I2.3ip obey the conditions needed in Lemma 

Proof. Unitarity of u and its periodicity in k 2 are obvious, while the CS symmetry is a consequence 
of Lemma 12.81 Also, (I2.29P is a direct consequence of (12.351) . 

Now let us consider the vectors in (12.311) . The only thing left to be proved is that they are the 
same when x = ±1/2. We have: 


N 

^'m(x, k2) = ^[^(a:, k2)~^]sm'^s{x, k2). (2.36) 

S = 1 


Then: 


N 


N 


N 


,(1/2, k2) = ^[m( 1/2, fc2)]rj^'r(l/2, fca) = ^[u(l/2, k2)\rj ^ ,9^,(^ 2 )^„^(-1/2, fca) 
r—1 r—1 m—1 

N N N 

= ^[u(l/2, k2)]rj ^ fdmr{k2) ^[u(-l/2, fc2)“^]smSs (-1/2, ^ 2 ) 


r—1 

N 


m—1 s=l 

= ^[u(-l/2, fc2)“^,5(A:2)u(l/2, fc2)]sjSs(-l/2, ^ 2 ) 

= S,(-l/ 2 ,fc 2 ), 

where in the last line we used (12.291) . 


(2.37) 

□ 


The last difficult technical problem we still have to solve is the construction of the approx¬ 
imation given in (12.321) . When this is done, (I2.34|) provides the unitary we are looking for in 
Proposition 12. 16 I which finishes the construction. 

The following lemma solves the approximation problem and is of independent interest: 
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Lemma 2.18. Let f3{k) be a continuous family of unitary matrices, Tj-periodic and satisfying the 
CS' property. Then there exists a sequence of real analytic families /3„(fc) of unitary matrices, 
'L-periodic, with nondegenerate spectrum at 0 and 1/2 and satisfying the CS' symmetry 
such that: 


lim sup ||/3„(/c) - /3(fc)|| = 0. (2.38) 

n^oo 

Proof. Before giving the proof, let us remark that each satisfies the conditions from Subsection 
12.7.41 hence they can be written as Pnik) = where each hn is continuous, periodic and with 

CS symmetry. This implies (I2.32|) since ||/3n(fc) —/3(fc)|| = ||Id — 

Now let us prove the lemma. We split the proof into two parts. First, we construct a continuous 
family /?„, with a completely nondegenerate spectrum at 0 and 1 / 2 , Z-periodic and with CS' 
symmetry. Second, we show how it can be made real analytic by preserving all other properties. 
Step 1. Assume that the spectrum of /3(0) consists of 1 < po < fV clusters of degenerate 
eigenvalues labeled as {Ai(0),..., App(O)} in the increasing order of their principal arguments. If 
s > 0 is sufficiently small and \k\ < s, due to the continuity of /3(-) we know that the spectrum 
of /3{k) will also consist of well separated clusters of eigenvalues. Let Cj be a simple, positively 
oriented contour, independent of |A:| < s, enclosing the j’th cluster and no other eigenvalues. Let 
llj{k) be the spectral projection of /3(x) corresponding to the j’th cluster. We have: 


n,(fc) = — / im - z)-^dz, ‘n,(fc) = n,(-fc). 


27r 


tCi 


lim IlnAfc) — n, 

fc^O 


= 0 , 


|fc| < s. 
(2.39) 


The matrix f3(k) is block diagonal with respect to the decomposition nj(fc)C^, i.e. 

/3(fc) = nj (fc)/3(/c)n, (A:). Define 

Po 

P(k) := £A,(0)n,(fc) = e*5:,"aiArg(A,(0))n,(fc)^ |^| ^ 

i=i 

Clearly, (3{k) is unitary, commutes with j3{k) and ‘/3(fc) = f3{—k) if |fc| < s. Define 'y{k) := 
fi~^{k)(3{k); we have that 'y(k) is unitary, commutes with /3(fc), *7(fc) = li—k) if |fc| < s and 
limfc_>o 7 (fc) = Id. Using the Cayley transform as in Subsection 12.7.31 see formulas l|2.26l) - (|2.27|) 
with (f>o = '^1 we can find a self-adjoint matrix hifk) such that 

7 (fc) = h{0) = 0, *h{k) = h{-k), [Uj{k),h{k)] = 0, \k\ < s. (2.40) 

Putting everything together and using that /3(fc) = fi{k + m) for all m S Z we obtain: 

I3{k + m) = Arg(A,(o))n,(fc)-eli(fc))^ ^ ^ (2.41) 

We proceed similarly for k € {|A: ± 1/2| < s}. Notice that due to the periodicity of j3, there will 
be the same number of eigenvalue clusters at k = ±1/2 and we can choose the same ordering for 
them, i.e. \j{\/2) = Aj(—1/2) for all 1 < j < Pi /2 < fV; note that there is no connection with 
the numbering at fc = 0. From the Riesz integral formula, choosing a contour Cj around each 
\j[\/2) = Aj(—1/2) and using that j3{k) = j3{k + l) we obtain nj(—1/2) = nj(l/2) and moreover: 

%(k) = n,(-k), %(l/2)=n,(-l/2) = n,(l/2), , |fc±l/2 |<s. (2.42) 

As in (12.411) we get: 

/3(fc±m) |fc±l/2|<s, Vm € Z. (2.43) 

In particular, nj(±l/2) are real and symmetric matrices. 
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The main idea is to locally perturb P around all integers and half-integers such that the 
spectrum is nondegenerate at these points and all symmetry properties are left unchanged. Let 
us now describe the perturbation. First, consider k = 0. Since 11^(0)* = 11^(0) = *nj(0), 
these matrices are also real and symmetric. Accordingly, Ran(nj(0)) is invariant under complex 
conjugation and we can find a real orthonormal basis. Let 1 < rj{0) < N he the multiplicity of 
Aj(0), i.e. the dimension of Ran(nj(0)). We have: 


n,(o) 


rj (0) 




dim Ra.n{Pj^i. (0)) = 1, 


Pj,L=P*>- = 


n 1 - 
Ji^3 


t p. — p2 


Define: 


(2.44) 


rj(0) 

MO) ■■= E - 1) (0), 1 < J < Po- (2.45) 

lj—1 

By construction, we have *Aj(0) = Aj(0). Seen as an operator acting on Ran(nj(0)), the spectrum 
of Aj(0) is nondegenerate and equals {0,1, ...,rj(0) — 1}. 

Now consider k = ±1/2. Since nj(l/2) = nj(—1/2) are real for all 1 < j < Pi /2 (see (12.421) 1. 
we can mimic the construction in (12.451) and obtain: 


A, (1/2) 


'■Al/ 2 ) 

E M - 1) Pj,m^ (1/2), 1 < J < Pl/2- 


(2.46) 


Now let ( 7 s : R !-)• [0,1], continuous, even {gs{k) = gs{—k)), (7s(0) = 1, supp((7s) C [—s/2, s/2]. 
Define: 


( Po \ /Pl /2 

E^ijXo) ±sE gs{k + m + 1/2) 
j—1 J ttiGZ yi—1 

From (12.411) and (12.431) we see that P{k) = on the set: 

Ds := {\k — mj < s} U {j/c — m — 1/2| < s}. 

m^Z 


(2.47) 


We see that Dg is symmetric with respect to the origin and consists of a union of disjoint open 
intervals centered around all integers and half-integers. Also, the support of Vs is contained in 
Dg, it is periodic and has the CS symmetry. Define /3g(A:) to be if fc g Dg, and let 

/3s (fc) = P(k) outside Dg. The family /3g obeys all the required properties and converges uniformly 
in norm to /3 when s —^ 0. 

Step 2. Here we will show how we can make /3g real analytic without affecting the other required 
properties. In order to simplify notation, we drop the subscript s and we assume that /3(.) is 
continuous, with CS' symmetry, Z-periodic and with completely nondegenerate spectrum at 0 
and 1/2 (hence at all integers and half-integers). Let f{x) = 7r“^(l ± If z/ > 0 we define 

fv{x) := v~^f(x/v). The function / is positive, symmetric f(x) = f{—x), f^, is an approximation 
of the Dirac distribution and is analytic on the strip {jlm(fc)j < v}. Then the matrix family 


Mk) := [ fAk - k')P{k')dk' 

3r 


(2.48) 


is real analytic and Z-periodic. A-priori jivik) is neither self-adjoint nor unitary, but it obeys 
*g,i,{k) = g,^(—k). Also, Pi/(-) converges uniformly to /3(.) when —)> 0, which implies: 

lim sup \\fii,{k)fi,^{k)* — Idjj = 0. (2.49) 

^->■0 feeR 
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If V is small enough, the operator is self-adjoint and positive, hence we can define: 

/3.(fc) := {^^,{k)^l,{kr}-^^^^i,{k), (]Ak)P.{kr = Id. (2.50) 

Clearly, /3[/(') is unitary, Z-periodic and converges uniformly to /3(-) when i/ —^ 0. This guarantees 
that f3i, has completely nondegenerate spectrum at 0 and ±1/2 for z/ small enough. We still need 
to prove that has the CS' symmetry and it is real analytic. 

If ly is small enough we can write: 

I3iy{k) = -^ [ z~^/‘^{^iy{k)^i,{k)* - z)~^^iu{k)dz (2.51) 

27’" J|z-l| = l/2 

Now the matrix family 

7.(A:) := /r,(A:)/r,(fc)* = [ U{k - k')f,{k - k")P{k')p{k"ydk'dk" 

has a holomorphic extension to the strip {|Im(fc)| < v}. If v is small enough, due to (I2.49P 
it follows that {fi^{k)fi^{k)* — z)~^ exists for all 1^ — 1| = 1/2 and is uniformly bounded on 
{|Im(fc)| < a} with 0 < a -C 17. Together with (I2.48|) this proves that has a holomorphic 
extension to {|Im(fc)| < a}. 

The last thing we need to prove is ^fdyk) = I3v(—k), which would imply the CS' property. 
From (I2.50p we see that if v is small enough one can write: 

l^vik) = ^ ar,[iiv{k)^yk)* - Id]"'/r^(fc), (1 ± ^ 

n>0 n>0 

We have already seen that ''^v{k) = /r^(—fc) and this leads to *'[^J-v{k)*\ = fi^{—k)*. Then: 

*{[fi^{k)n„{k)* - Id]”/Xj,(fc)} = ^i„{-k)[fj,^{-k)*n^{-k) - Id]” = [fi^{-k)n^{-k)* - Id]”/i^(-/c) 

where the second identity can be proved by induction with respect to n. The proof of the lemma 
is over. 

□ 


3 Proof of Theorem 11.11 

We will use the generic notation a and C for a finite number of positive constants which will 
appear during the proof. 

The difficulty stems from the fact that since Pb is not norm continuous in b, one cannot use the 
continuity approach of [NN2] . The first step in overcoming this difficulty is to use an old ansatz by 
Peierls [Pej and Luttinger [Lu] and to couple it with the regularized magnetic perturbation theory 
in order to construct (out of the Wj^-y’s) an intermediate orthogonal projection Ilf, such that 

lim||nf,-Pf,|| = 0. (3.1) 

h—>0 

The second step is to construct a localized orthonormal basis for Ran(nf,), while the last step is 
to transfer it to Ran(Pf,) by the ‘continuity’ method of [N21INN2] . 

If |5| is small enough, the orthogonal projection Pb can be written as a Riesz integral (see 

Pb = ^ f {Hb - z)~^dz 

J 0 

where C is a simple, positively oriented contour which surrounds do (hence Ub) and no other part 
of the spectrum of Hq. By writing 

Pb = ^ j [{Hb - z)~'^ - {Hb - i)~ydz = ^ J {z - i){Hb - z)~^{Hb - i)~^dz 
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and using (I1.13L one can prove ( see also Section 4.1 and 4.2 in [C]) that Pi, has a jointly 
continuous integral kernel Pt,(x, x') and there exist a > 0 and C < oo such that 

|P6(x,x')| < Vx,x' e |6| < bo. (3.2) 

In particular, this holds true for 6 = 0. Since Po(x, x')wj_.y(x')dx', we have that 

each Wj^.y is continuous and by applying the Cauchy-Schwarz inequality and (13.21) we obtain the 
existence of some constants C < oo and a' < a such that 

sup e“ < C'. (3.3) 

xeR2,7Gr 


This implies the identity: 


N 

Po(x,x') Wy^(x)w—(x') 

3=1 7er 


(3.4) 


where the above series is absolutely convergent; let us show why. The estimate (j3.3|) can be 
reinterpreted as |wj,.7(x)| < for some C < oo and a > 0, uniformly in 7 and j. Using 

the triangle inequality we obtain: 


N 

EE |i(;j_.y(x)| |wj_.y(x')| < const e ’^11/2^ 

3=1 76r 


(3.5) 


Later we will need the following bound which holds true due to our sparsity assumption on T: 

sup ^ ||c — < 00, Vfc > 0. (3.6) 

cgr=:^ 


3.1 Proof of (i). 

Applying the Riesz integral formula in (I1.13|) we obtain: 

Pb(x,x') - e*'’-^(’^’’^')Po(x,x') < C |6| e 


^ —O' X —X 


(3.7) 


Define: 

P,(x,x') :=P''^(’''’'')Po(x,x'), (3.8) 

and denote by Ph the corresponding operator. Applying a Schur-Holmgren estimate using (IMI) 
leads us to: 

||Pb - P&ll < U |6|, where |Pb(x,x') - P6(x,x')| < C |6| H. (3.9) 


Define the kernel: 

N 

Pb{i^, x') = ^ ^ e**‘^^^’’*'^Wj^.y(x)e®*"^(’‘'’'>')wj, 7 (x'), (3.10) 

3=1 7er 


and by Pb the corresponding operator. The series is absolutely convergent (see the estimate in 
(13.51) 1 and defines a jointly continuous function. Moreover, we have: 

|Pb(x,x')| < (3.11) 

From the definition of (()(x, x') (see (11.141) and (11.111) 1 it follows that it equals the magnetic flux 
through the triangle with corners situated at 0, x and x'. 
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Let /Z(x, y,x') denote the magnetic flux through the triangle with corners situated at x, y 
and x'. Using Stokes’ theorem and the fact that |i?(x)| < 1 we have: 

/l(x,y,x') = (j){x,y) + cj){y,x.') - (j){x,x'), |//(x,y,x')| < ^|x-y| |y-x'|. (3.12) 

Using the various definitions, the triangle and Cauchy-Schwarz inequality we obtain: 

N 

\Pb{x,x') - Pb{x,x')\ < y]y] 2 |sin( 6 /Z(x, 7 ,x')/ 2 )| |wy 7 (x)| |wj-^(x')| 
j=i 7 Gr 

< sup V ||c - 

which leads to the existence of a C < oo and a > 0 such that 

|Pb(x,x')-n(x,x')| \\Pb-Pb\\<C\b\. (3.13) 

Consider the closed subspace 

Sb ■= closure{Span{e*^'^^'’'’'^ri;j^.y(-) : 7 S L, j e {1, 2,..., N}}} 
and its associated Gram-Schmidt matrix 

Mb{j,j;i',f) ■■= ■ (3-14) 

Lemma 3.1. Consider the self-adjoint operator Mb acting on the space /^(L) (8> given hy the 
matrix elements Mb{'^,j',^',j'). If bo is small enough, then there exists some a > 0 and C < 00 
such that if \b\ < bo: 

V - Mbi^,r,y,f)\ < C \b\ e-“ll^-'^'ll. (3.15) 

Moreover, Mb > 1/2, and if (3 & {1/2,1} then 

- M-^(7,j;7',/)| < C \b\ ( 3 . 16 ) 

Finally, the vectors: 

Y. E[M-^/^](7',/;7,J>*^^("’^')^y.7'(x) (3.17) 

7'ez2y=i 


form an orthonormal basis of Sb, being uniformly exponentially localized around 7 , i.e. 

|^'j. 7 ,b(x)| < C |^'j, 7 ,b(x) - e*'’‘^(^’'^)w;j,.^(x)| < C | 6 | (3.18) 

Proof. Define Db := Mb — Id. Since Sjj>6-fj' = {wj','y'i')\wj^-f{j )and = 0 we have: 

Dbh,j;i',j') = [ Wj',7'(x) Wj,7(x) 

JS.'i ^ •' 

and 

|Pb(7,j;7',j')l < [ ||7'-x||e“'ll'^'-"ll|a;,7y(x)| ||7-x|| |«;..7WMx 


for some small enough a' > 0. Using the Cauchy-Schwarz inequality we find: 

|P6(7,J;7,/)I < \b\ 


(3.19) 
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If n > 2 we have 


jl,7l jn-l,7n-l 


and (j3.19l) implies: 


Vfc.r,7"^ / 

< l&rC”, Vn > 1. (3.20) 

The estimate (13.191) also implies that Db goes to zero with b in the operator norm, hence Mb is 
close to the identity operator, thus M^^ = (Id + Db)~^ exists and can be expressed as a norm- 
convergent power series around zero. Then (13.201) implies (13.161) . Finally, it is straightforward 
to check that the vectors defined in (13.171) are orthogonal and span Sb, while (13.181) is an easy 
consequence of (13.161) and the triangle inequality. □ 


The orthogonal projection associated to the subspace Sb is 

N 

(3-21) 

7eri=i 

Under the conditions of Lemma [O] using (13.171) and the self-adjointness of Mb we obtain that lib 
has a jointly continuous integral kernel given by: 


n5(x,x') = 


N 

E E 

7',7"Gr j'= l 




where the series converges absolutely. Using p.l6p and (I3.10p we find that there exists some a > 0 
such that 

|nb(x,x') -Pb(x,x')| < C |6| (3.22) 

Using (13.9p and p.l3l) we obtain: 

|nb(x,x')-Pb(x,x')| ||nb-Pb|| ^Cl^l, (3.23) 


which finishes the proof of (13.11) . 

The next step is to consider the Sz.-Nagy unitary operator Ub which intertwines Pb with lit,, 
i.e. PbUb = UbUb, given (as in (I2.4|) l by the formula: 

Ub := {PbUb + (Id - P6)(Id - nb)}{Id - {Ub - Pb)2}-^/^ (3.24) 


A consequence of (I3.23P is that for small enough &, the operator Ub — Id has a jointly continuous 
integral kernel which obeys the same estimate as in (13.231) . Then the set of vectors: 


N 

:= Ub'fj,,,b = E E (3.25) 

7'er j'=i 

form an orthonormal basis of Ran(Pb) satisfying (11.151) and (11.161) . 
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3.2 Proof of (ii). 

Remember that here the magnetic field is constant with magnitude -B(x) = 1 and F is a periodic 
lattice embedded in In this case, the phase is not just antisymmetric but it is also bilinear, 
since now ^(x, x') = ^(x{x 2 — X 1 X 2 ) or, with the usual abuse of language, (p(x, x') = • (x' Ax). 

The identity (x — y) A (x' — y) = —x' Ax + yAx + x'Ay implies: 

</)(x, y) + x) - x') = <^(x' - y, x - y) = <^(x' - y, x - x). (3.26) 

Denote by the unitary magnetic translation given by: 

[T,V](x) = e*'"^(-’^)/(x-7). 

Using (I1.17L p.26l) and (13.141) we obtain: 

= (7-y(^«A.o)|T^(u'yo))L2(R2) =: Mb(rt - i] j, j'), 

where 

Mb{l"]j,j') = f ^Wj,o(y-7")''«A.o(y)dy = (u;jyo|Ty,('u;j,o))L2(R2). (3.27) 

JR2 


_2 /2 

Lemma 3.2. There exists some bo small enough such that the operator ' has a matrix of 
the form 

,j') -■ e''’‘^^^’^^7;(7-7';j,/), 
where \Tb{y"]i,j')\ < C'e““l'^ I for every 0 < b < bo- 

Proof. In the proof of Lemma dSH) we introduced the operator Db = Mb — Id in /^(F) (g) and 
showed that if |&| is small enough then = (Id + D^;)”^/^ exists as a power series in Db- The 

matrix elements of Db inherit the symmetry properties of My. 


'^bh";jj') := Mbij'-Jj') - Vo%A Db{'^,j\i,j') = )t>&( 7 - 7';j,/). 


Then using (|3.26p and the fact that T is a lattice we obtain: 


N 

j"=i 7"er 

= e*f-b(7.7') ^ ^ 
j"=i7"6r 

=:e*''^(7-7')Df (7-7';j,/), (3.28) 


N 


where 

b'\ — JHiin")'! 

j"=i7"Gr 

If n > 3 one can prove by induction that ^^( 7 , j; 7 ', j') := where: 


(7;a/)= E E -i"-,j.]")Vb(,i'-,j".]’). 




N 

E E 

"=i7"Gr 






This structure will be preserved for any convergent series in Db, while the exponential localization 
is a direct consequence of (I3.16F 


□ 
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(3.29) 


Using the above expression for j;7^ jO (13-171) we obtain: 


7'er j'=i 


*bb(7'.7)7-^(y _ _ y). 


We note the identity 

<(’(7 - 7', X - 7') = <(’(x - 7 ,7' - 7) 
which together with (13.2611 gives: 

</>(x, 7') + <('(7', 7) - 7) = 0(7 - 7', X - 7') = «(>(x - 7 ,7' - 7)- 

Replacing this in (|3.29p we obtain: 

N 

d/,.^,,(x) = ^ ^ ryy - 7;/, - 7') 

7'er 

N 

= E E ry7";/, - 7 - 7") 

7"6ri'=i 


AT 


Denote by: 


7 "er j'^i 


N 

'^0,b-= E Ty'iwfA- 

7"eri'=i 


(x). 


(3.30) 


(3.31) 


From (13.301) we have = T^{uj,b)- The last technical result we need is the following: 

Lemma 3.3. Both Pb and Ilh commute with the magnetie translations and so does Ub- 

Proof. Since Hb commutes with y, so does Pb- We only need to prove this for lib. Using the 
formula = T^(uj,b), we can express the kernel of Hb (see (13.211) as: 


N 


nb(x,x') = E E® 

7'er j=i 


z60(x,7')y^. b(x - - 7')- 


We have: 


N 


[T^^nb](x,x') = ^ - 7 )e-*''^("''^'y,-b(x' - 7 ') (3.32) 

7'eri=i 


and 


N 


[nbr^](x,x') = J2 Xle*'’‘^^''’^'^Wj-b(x-7')e-*'"^(^'+T''T'')e*'’'^(’''+T'-T'yj-b(x'-7' + 7). 

7'Gr j=i 

Changing 7 ' with 7 ^ + 7 in the second formula we obtain: 

N 

[nbT^](x,x') = E - 7 ' - b(x' - 7 '). (3.33) 

7 'er j=i 

Using the bi-linearity and antisymmetry of (j) we see that the magnetic phases in (I3.32[) and (I3.33|l 
are the same and the two kernels coincide. □ 
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Denote by Uj^b ■= UbUj^t- Since commutes with Ub, we obtain that 

— Ub^j^j^b — UbT^i^Uj^b') — 

and: 

N 

Pb = \'’''y{'^0,b)){T^{Uj^b)\^ ^j,0,b = Uj^b- (3.34) 

j=i 7Gr 

In particular, this proves (11.181) . The last thing we need to prove is (11.191) . or equivalently: 

Uj^b = Uj-b- (3.35) 

This last identity would be implied by two others: 

uj^ = uj-b and C7^= C/-b7, V/e L^(R^). (3.36) 

In order to prove uJJ = Uj^-b, we see from (I3.31|) and (11.1711 that the only thing we need for this 
is 

WWJ) = T-b{l"-,3',3). (3.37) 

From (j3.27l) and (11.1711 we obtain that ■Mb{Y';j',j) = A4-b{'y";j',j), hence M^i/; = M-bij^ for all 
■0 G P{T)®<C^. By functional calculus, the same identity is obeyed by which proves (13.371) . 

hence uJJ = Uj^-b- Since '^j,-y,b = ’T^i.Uj^b) we immediately get '^j,-y,b = '^j,-y,-b- From (13.211) we 
obtain that Ilf,/ = n_f,/ for all / G T^(R^). Moreover, since Hbg = H-b 9 for all g G (^((“(K^), a 
standard argument shows that Pbf = P-bf for all / G T^(]R^). Now using the formula (I3.24p we 
obtain the second identity in (|3.36l) and the proof of Theorem II.Ill is over. 
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